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1. Introduction 

To reproduce the physical space- space time out of some more fundamental variables, 
rather than introduce it explicitly, has been always a tempting idea in the quantum field 
theory. While the string theory sets up a fruitful framework in which the space-time is 
dynamically created out of the fluctuating coordinates of the strings, potentially not less 
fruitful may be the attempts to encode the space-time, together with the quantum fields 
themselves, into the dynamics of large fluctuating matrices. A quantum field theory or a 
string theory should be described in this case in terms of specific matrix integrals containing 
finite amount of matrices of infinite size. The earliest proposal of this kind belongs to 
T. Eguchi and H. Kawai |JJ, followed by a few important precisions and modifications 
' ' ' • I n these works the QFT with N x N matrix valued fields can be reduced in 
the large N (planar, or 't Hooft) limit to finite dimensional matrix integrals in the same 
limit. 

Another successful enterprise of this kind was a formulation and solution of non-critical 
string theories (associated with the two dimensional quantum gravity in the presence of 
some matter fields) in terms of U(N) invariant matrix integrals ||, 0. 

More recently, one of the most fruitful ideas of this kind in the superstring theory was 
the proposal of E. Witten || generalizing an old idea of Chan-Paton to describe the low 
energy physics of branes by various reductions of super Yang-Mills theory. The diagonal 
components of the effective SYM vector potentials play the role of space-time coordinates 
of branes there. 

The concept of the physical space-time built out of discrete degrees of freedom has 
become especially inspiring due to fundamental questions in quantum gravity, such as 
microscopic explanation of the thermodynamics of black holes. 

Our discussion in this paper will be confined to the matrix model formulation of 
quantum field theory. To set up the problem let us recall that the old Eguchi-Kawai (EK) 
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reduction reproduces correctly only the planar sector of a matrix field theory, whereas the 
non-planar corrections (1/N expansion) were never incorporated into this scheme, let alone 
a nonperturbative formulation of a QFT with finite number of field components (finite N) 
in terms of some matrix model of a matrix or matrices of infinite size. The reason for this 
difficulty is mostly due to the lack of reduced momenta running around nontrivial cycles 
of non-planar graphs in the EK reduction scheme. 

The aim of the present paper is to propose a new formulation of a finite component 
scalar field theory (finite N) in terms of an n x n matrix integral in the limit n —>■ oo. One 
can say that it incorporates all orders of the 1/N expansion and is in principle a possible 
non-perturbative definition of the original scalar field theory. For example, the usual scalar 
(f> theory can be formulated as a one matrix integral in external matrix sources. We will 
also show how to incorporate fermions into this scheme. Unfortunately, we did not find so 
far any natural way to formulate the four dimensional QCD in this way. 

Our construction is in some sense T-dual to the old EK scheme: we use the diagonal 
matrix sources of quenched coordinates instead quenched momenta of the EK scheme. As 
a consequence of it the original scalar field "lives" on the graphs dual to feynman graphs of 
our matrix model. To control the parameters of these graphs (say, to make them exactly 
(f) A graphs) we apply the methods worked out in and |TD| for the so called model of 
dually weighted graphs (DWG). 

We will be able to generalize our method to fermions and to their yukawa interactions 
with the scalars, but for the moment we don't know a natural way to introduce gauge 
symmetry into our approach. So to formulate QCD is an interesting challenge in our 
framework. 



1 Using this abbreviation we remember of course about the important modification of the 
original, not quite working, reduction scheme of Eguchi and Kawai related to the quenching of 
momenta Q and twisting Q 



2. Difficulty with 1/N corrections in the EK reduction scheme 

Let us remind the essence of the old (and unsolved, to our knowledge) problem of 
1/N corrections to the reduced version of planar field theory. We will mostly discuss a 
matrix version of scalar field theory in the euclidean four dimensional space described by 
the action 

S = N J d^xtx {{d^) 2 + V {<!>)) (2.1) 

where (f) n m(%) is an N x N hermitian matrix field and V(0) = X]fc>2 ^ k( f )k * s a sca l ar 
potential. 

The EK reduction (in the most natural, Parisi formulation [EI) goes as follows: take 
the following particular dependence of <p on the coordinates: 

4>mn(x) = e~^ Pm ' x "> (j) mn e^ Pn ' x \ no summation over m,n (2.2) 

The action then takes the form 

S = Vtr([p^ ( j ) } 2 + V( ( p)) (2.3) 

where V is the 4D volume of the physical space, p^ = diag(p^\ ■ ■ ■ ) are D diagonal 
matrices of quenched momenta scattered uniformly |4[] in a large 4D "momentum box" of 
a size A 4 , where A is UV cutoff, and the original functional integral over scalar field is 
replaced by a single matrix integral over ^-independent (f) mn . The planar sectors (leading 
large N approximation) of the two models, the original matrix scalar field theory in 4D 
and the reduced OD one matrix model, coincide. It is immediately clear from the double 
line representation of planar graphs in the reduced theory (fig.l): each face of such a graph 
is associated with a closed index loop and with the momentum variable pj/ carrying the 
same index i: each double line propagator Df-, = -L ^"'J^U depends on the difference of 
the momenta of adjacent faces. 
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Fig.l: A fragment of a planar diagram in the double line notation and the propagator 
depending on the quenched momenta pj-, Pj, k = 1,---,N in Eguchi-Kawai reduction 
scheme. 

The planar part of the free energy of the reduced model F p i an = UmN^oo^2 ^ogZ^ ( 
Zn is the partition function) can be schematically written in the following way: 



where J2 G goes over planar feynman graphs G of the scalar field theory, v,v we label 
respectively original and dual vertices of a graph G, i$l,---,N is the index associated 
with the dual vertex v and n<^«'> S oes over ^ ne edges < vv' > of the graph G dual to G. 
In the large N limit only planar graphs survive in both models and the sums over 



the planar sectors of the original matrix scalar field theory and of its reduced version are 
equivalent. This equivalence extends to any one point U(N) invariant physical quantities 
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(2.4) 



indices reproduce the integrals over 4D momenta: 



N 



I\p\<A d4 P---- HenCe 



of the type Ok =< jj(f) k (x) > but is known to fail for the multi-point correlators since 
< jj(p h (x)jj(j) l (y) >=< jj(p k (x) >< jj4> l (y) > plus 0(1/N 2 ) corrections, different in 
two models. The EK reduction fails to describe correctly the higher 1/N (non-planar) 
corrections to the original scalar field theory 

The reason for this failure is well known: we cannot represent all momenta running 
through the propagators as differences of momenta of the adjacent loops (faces) on the 
graphs of a non-spherical topology. If we did so (and it is precisely the case of the topo- 
logical expansion in the EK reduced model) the momenta running along topologically 
nontrivial cycles of a non-planar feynman graph would be zero. To see it one takes any 
nontrivial closed path on a dual graph (connecting dual vertices or original faces) and 
calculates the total momentum running through the propagators crossed by this path as 
(Pi ~Pj) + (Pj ~ Pk) + ■ 4 ■ + (pi — Pi) =0. For example, for the torus topology we have two 
momenta of the original 4D theory missing in the reduced version: they flow through two 
nontrivial cycles of the torus. 

The reason for the failure is simple but the remedy is not easy to find, at least in case 
of the EK reduction involving reduced momenta. 

In the next section we will show that the goal of construction of a matrix model 
describing a finite N (including the most frequent N = 1 case) scalar field theory can be 
achieved by introducing quenched coordinates instead of quenched momenta. 

3. Matrix model formulation of finite iV scalar field theory 

Now we will propose a one matrix model in external matrix fields which will be 
equivalent, at least perturbativly, graph by graph of any topology, to the original 4D 
finite iV matrix scalar field theory ( |2.1p . We will show that the free energy and, with an 
appropriate definition, the physical quantities of the field theory ( |2.1| ) at finite N coincide 
with those of the matrix integral over a hermitian n x n matrix $ in the limit n — > oo: 




(3.1) 
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with the action: 

S = NTr M $] 2 + HIM ~ A*)) (3.2) 

Here the matrix $ lives in the n = p ■ q dimensional vector space Af which is a direct 
product Af = V x Q of vector spaces of smaller dimensions p and g, correspondingly. I n 
is the n x n unity matrix. Both dimensions p, q go to infinity as n — > oo; and A are 
external (fixed) diagonal matrices of the form 

X^ = x„x I Q , (3.3) 

where x^ = diag(xy , • • • , zjL ), 

A = / P xa (3.4) 

where a = diag(a\, ■ ■ -a q ). 

So the matrices x M on the one hand and the matrix a on the other hand live in 
orthogonal subspaces. 

The matrices with /j, = 1,2,3,4 will play the role of quenched coordinates: the 

(i) 

points with coordinates should be distributed uniformly in the physical space box of 

1/4 

a size L 4 , which is the size of our system. The ultraviolet cutoff is defined as A ~ . 
Obviously the thermodynamic limit of (infinite volume) corresponds to p — > oo with A 
fixed. 

The matrix a will encode the information about the couplings of scalar potential 



in the form: 



k 

k>2 



t k = ~ tr Q a k (3.6) 



where the trace trg goes only with respect to the vector space Q. Here t2 = m 2 corresponds 
to the mass squared. 
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The parametrization (|3.6| ) of the couplings reminds the so called Miwa variables widely 
used in the theory of r-functions of the hierarchies of integrable differential equations. 
It is also used in the representation of characters of the group GL(N) through Schur 
polynomials (see [|10| for the details). 

Obviously the last formula can be in general true only in the limit q — > oo. Note that 
N is kept finite fixed parameter here. 

/.From ( |3.5|) and (|3.6|) the potential can be also written in the form 



^) = -|Emi-m ( 3 - 7 ) 

3=1 

In the limit q — > oo we can parameterize in principle any potential (including a polynomial 
one) by such a sum of logarithmic terms, but a^s need not be necessary real. They can 
be taken, say, in complex conjugate pairs (the potential becomes even in this case). 

For instance we can choose in such a way that in the limit q — > oo they will 
reproduce any polynomial potential (|3.5| ): 

a,j = e i9j (3.8) 



where 



27r 2N ^ tm ■ Znjm 

Oj = — j H } — sm (3.9) 

q pq tu q 



Indeed, for k > we have 



2in 



E,,, I , 2A^ \ ^ t m . 2irjm \ N 
e lke > ~ p , q ^oc V e 5 Jfc 1 + ifc V — sin = t fc . 

j=l j=l \ FH m>l H ' 



The proof of the equivalence of the QFT ( |2.1| ) and the matrix model ( |3. 2|) is very 
simple. Let us consider any feynman graph of the theory ( |3.2|) (dotted line on fig. 2) 
together with its dual graph (solid line on fig. 2). The matrix structure of the theory 
prescribes to use the double line notations for the propagators, so such graph has a fixed 



topology with the genus defined by the Euler formula. Each single line carries now a double 
index corresponding to the product of spaces V x Q. 




Fig. 2. Original and dual graphs and the dual propagator depending on the quenched 
coordinates xj-,Xj in our one matrix integral representation of scalar QFT. We show the 
pairs of indices 6, k and c, j belonging to Q, V spaces, correspondingly, running around 
original faces, and hence placed at dual vertices. 

The propagators of the original graph are given by: 

1 a b a c 



N (xi - Xj)' 



:0~bb'0~cc'0~ii'0~j 



Here we attributed the factors to the propagators rather than to the vertices which 
can be achieved by the change of the matrix variable: $ — > A~ 1 / 2 QA~ 1 / 2 . It is easy 



to see that we obtain at each face of the original graph a weight tiQa k where k is the 
order of this face (number of edges). It happens in the same way as in the so called 
matrix model of dually weighted graphs (DWG) - a one matrix model with the action 
S = ntr[$ 2 + W(A$)] studied in and solved in [fToH : apart from original couplings 



coming from the potential W(z) we also obtain in the DWG model the dual couplings tj~ 
weighting the faces of different orders k (or, which is the same, the vertices of the dual 
graph with the coordination number k). 

The factors ^ x .l x .y can be now attributed to the dual propagators (crossing the 
original ones) and the indices i,j are running around the faces adjacent to the original 
propagators or, which is the same, attributed to the corresponding dual vertices (see fig. 2). 

It is natural to formulate the result for the free energy in terms of feynman expansion 
with respect to the dual graphs G. It can be written in the following way: 

^E^nt^t) n K'-^r 2 (3-io) 

Q V lE=l <V'S"> 

where goes over all dual graphs of the matrix model ( |3.2|) , g is the genus of a graph G, 
Y[y goes over all vertices v of this graph with coordination numbers k, tk = — jy Y^j=i a j 
are the couplings attached to these vertices and n<^'^»> S oes with respect to all edges 
of G connecting the vertices v and v' . At each dual vertex v there is a sum taken with 
respect to the index i% corresponding to the subspace V . 

Note also that due to the specific logarithmic form of the interactions in ( |3.2j ) 
-Trln (I-A&) = J2 k ^{A®)* all faces of dual graphs appear weighted with the factor 1 
(the factor v compensates the cyclic symmetry of each dual face). So we see that (|3.1U| ) is 
given by the sum over connected graphs G of all genera waited by N 2 ~ 2g , with unrestricted 
face order and with the vertices waited by couplings tk % = jr trg d fcc . Hence these are 
precisely the original graphs of the scalar matrix field theory ( |2.1|) with the appropriate 
4D massless propagators in the coordinate space. The mass m is taken into account by 
the presence of the coupling t<i = — m? in the scalar potential. 



It is left to add that the summations over the indices i% = 1, • • • ,p can be substituted 
by the integrations in the large p limit: 



1 p r 

-p^x, / dV-. (3.11) 



y i=l 

In this way we can reproduce the correspondence between the matrix scalar 4D QFT 
( p.l|) and the zero dimensional matrix integral (|3.2|) , graph by graph of any topology. 
Hence they coincide, at least in any order of the perturbation theory. 

Let us stress again that N is a fixed parameter in our construction. It does not even 
need to be integer, although integer iV seems to be singled out by the fact that we can 
represent the corresponding determinant as an integral over an N- vector of complex n x n 
matrices Mi, I = 1, • • • , N: 

r N 

exp[-iVTrlog(i"- A$)] = J[d 2N2 Mi exp[-TrM+(J - A$)Mj] 

J i=i 

In particular, for N = 1 ( |3.2| ) is equivalent to the usual one component scalar field theory 
with the action S = J d A x[{d^f + V (</>)]. 

To calculate the one point correlators we use the following correspondence between 
the averages in the scalar QFT and its matrix model (MM) formulation (we take N = 1): 

< J d 4 X(j) k (x) >QFT= < tw(A$) fc >MM (3.12) 

For the two point correlator we have the correspondence: 

< <j>{x)(j>{y)) > QFT = J^T^i < (^(Atyh (tr jsf(A$)) j:j > M m (3.13) 

where the traces are taken only with respect to the subspace Q and the matrix indices ii 
and jj of the subspace V are fixed and chosen in such a way that (xi — Xj) 2 ~ (x — y) 2 . 
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4. Reduction in the presence of fermions and yukawa interactions 

Although we don't know any natural way to build a reduction of Yang-Mills theory 
for finite N, the scalar field theory ( |2.1|) is not the only interesting QFT which can be 
reduced to a matrix model in this way. 

Let us consider as an example of application of our method the reduction of the QFT 
of massless Dirac fermions and massless bosons with yukawa interaction in 4D with the 
action: 

S = Jd 4 x j(d^) 2 + ^ - (4-1) 

We stress that here <p(x) and i/j(x) are usual bosonic and dirac fields: has only one 
component and ip is a dirac spinor. The reduced version of this model is given in terms 
of matrix integral over 2 hermitian n x n matrices $ and \& in the auxiliary linear space 
Af = QxV 

Z = e F = [ cT 2 $cT 2 # e~ s (4.2) 



with the action 

S = T W (LY M , $] 2 + [X^ [X v , m 2 - tT D 1h(/at - 7m x [X» + Intfv - A$AW)) 

(4.3) 

Here tr^> is taken with respect to the dirac indices of 7^ matrix which is in the direct 
product with all other matrices in the third term. The A matrix is chosen in such a 
way that p^ =1 = \5k,3 for k > in the large q limit. The last can be achieved for 

/ n1/3 

example by the following choice of a: aj = f |^ J e l0j with 9j given by the equation 
= 6 + I sin 36*. To prove it it is enough to calculate the density of 6>'s in the large q 

limit: p{6) = ±§f = ^(1 + | cos 30) which gives tr Q a fc = j\e ike p(6) = A5 M for k > 0. 
All other definitions are the same as in the previous section. 

To verify the perturbative equivalence of the QFT ( [4.1| ) and the matrix integral (|4.2| )- 
( |4.3|) we compare again the feynman graphs of the former to the dual feynman graphs 
of the latter. Due to the choice of the matrix a dual graphs of ( f4.3| ) will contain only 
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triple interaction vertices and due to the log type potentials the order of dual faces will be 
unrestricted (see fig. 3). 




fermionic propagator 
bosonic propagator 

dual fermionic propagator 
dual bosonic propagator 



Fig. 3: Original and dual graphs of the QFT of scalars and fermions with yukawa 
interactions (in single line notations). 

Note that feynman graphs of the QFT ( |4.1[ ) have two types of faces: fermionic loops 
built out of Dirac propagators and the loops with the boundary built from interchanging 
dirac and bosonic propagators. The first type of loops will be generated by the third 
term in ( |4.3| ) and the second kind - by the last term in (|4.3| ). The yukawa interaction 
vertices correspond to the loops of dual graphs of the matrix model ( |4.3| ). Note the 
(— ) sign in front of the third term which gives the fermionic statistics ((— ) sign for each 
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fermionic loop). Each fermionic loop of k vertices is equipped, as it should be, by the factor 
trr>(7 Ml • • - 7/i fe )- Finally, the massless 4D fermionic propagator |^_^|4 will appear as a 

/ jJ, u \ 

function of difference of quenched coordinates f 2_ — -rf^. 

Hence in the same limit p, q — > oo as for the scalar QFT of the previous section we 
reproduce graph by graph (with the weights independent on the topology of graphs) the 
free energy of the QFT ( flip out of the matrix integral ( |4.3| ) . 



5. Conclusions and comments 

We proposed a reduction of 4D quantum field theories with a finite number of fields 
(finite N in case of matrix fields) to a matrix integral over infinite matrices. The phys- 
ical four dimansional coordinate space is encoded into the components of D=4 auxiliary 
diagonal matrices of quenched coordinates. The reduction is different from the old Eguchi- 
Kawai reduction with quenched momenta: it reproduces correctly not only the planar 
approximation but also, at least perturbatively graph by graph, all non planar corrections. 
The coordinate space is introduced in our reduction scheme in a way which reminds the 
description of coordinates of D-branes in H. The large n limit that we use is different 



from the usual 'tHooft limit and rather similar to the one adopted in |12[ for the matrix 
model of M-theory. 

We showed that fermions can be also naturally introduced into this reduction scheme. 
A few comments are in order: 

1. The QFT's in dimesions different from D=4 can be also reduced in this way 
but since the scalar propagator is different from l/(x — y) 2 the corresponding matrix 
representation of it looks "nonlocal" (i.e., different from Tr[AT M ,$] 2 ). One can use the 
formula TrfX^, [• • • , [X^ k , $]] • • -] 2 = ^\ -{xi — Xj) 2k \$ij\ 2 in the action of the reduced 
matrix model to supply the dual graphs by any propagator D(xi — Xj). 

2. Although our new matrix formulation of some QFT's hardly helps for solving them 
analytically it may provide a new numerical approach to their study: our matrix model 
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does not deal with any 4D lattice and the approach to the thermodynamical limit might be 
much faster then in the conventional Monte Carlo algorithms based on the lattices. One 
could also envisage some real space renormalization schemes where the renormalization 
flow would be considered with respect to the size of matrices. 

3. The authors of the papers ||13|| , JTJ] propose a formulation of the non-comutative 



QFT using reduction to the matrix models with quenched momentum matrices P^ obeying 
the Heisenberg commutation relations: [P^,P V ] = iB^ u where is an antisymmetric 
D x D matrix of C-numbers. For our model, a natural non-comutative generalization 
would occur if we take the same commutation relations for the coordinate matrices 
[Xfj,,X v ] = iC^ v . Since the coordinate and momentum matrices are connected in |13[ by 
the relation P^ = B^yXy it is natural to expect that we get the same non-comutative scalar 
field theory if we take the matrix C = B . We haven't yet a proof of this statement. 

Let us also mention the most obvious problems and questions concerning our matrix 
reduction of QFT's: 

1. We did not manage to find a natural matrix reduction of QCD with 3 colors. In 
principle we could try to do it in a particular gauge by systematically reproducing all 
feynman diagrams by the method we presented in this paper. But it would be much more 
instructive to find some general principle for such reduction arising explicitely from the 
gauge invariance. For the moment such a principle is missing. 

2. It would be interesting to find a modification of our approach similar to the twisted 
EK model proposed in ||. In || the momentum space appears as a classical vacuum of a 
large N lattice QFT with modified couplings. We can imagine a similar twisting of reduced 
matrix model with the classical vacuum solution generating the coordinate space. 

3. A related but more ambitious question: can we build realistic models of fundamen- 
tal interactions as some reduced matrix models where the physical space would emerge 
due to some symmetry breaking procedure? In other word, can our World be described 
by a Matrix Model? We know some of the attempts of this kind in the superstring theory 
and M-theory [12], [15[j. Our construction might be useful to approach this problem from 



a different direction. 
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